Nutrient Shielding in Clusters of Cells 
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Cellular nutrient consumption is influenced by both the nutrient uptake kinetics of an individual 
cell and the cells' spatial arrangement. Large cell clusters or colonies have inhibited growth at the 
cluster's center due to the shielding of nutrients by the cells closer to the surface. We develop an 
effective medium theory that predicts a thickness £ of the outer shell of cells in the cluster that 
receives enough nutrient to grow. The cells are treated as partially absorbing identical spherical nu- 
trient sinks, and we identify a dimensionless parameter v that characterizes the absorption strength 
of each cell. The parameter v can vary over many orders of magnitude between different cell types, 
ranging from bacteria and yeast to human tissue. The thickness i decreases with increasing v, in- 
creasing cell volume fraction (j), and decreasing ambient nutrient concentration tpao- The theoretical 
results are compared with numerical simulations and experiments. In the latter studies, colonies of 
budding yeast, Saccharomyces cerevisiae, are grown on glucose media and imaged under a confo- 
cal microscope. We measure the growth inside the colonies via a fluorescent protein reporter and 
compare the experimental and theoretical results for the thickness i. 
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I. INTRODUCTION 

Nutrient uptake is essential for all life and has been 
studied in a variety of model organisms; however, the 
physical mechanisms involved in the uptake are not yet 
well understood. A cell commonly takes up nutrients 
from its surrounding medium via facilitated diffusion or 
active transport: Specialized transporters on the cell sur- 
face move the nutrient down (facilitated diffusion) or up 
(active transport) a concentration gradient from the am- 
bient medium into the cell. The cell's absorption of nu- 
trients in its immediate vicinity sets up and maintains a 
concentration gradient outside the cell that gradually de- 
pletes the nutrient from the surrounding medium. This 
is the method used by budding yeast cells for acquiring 
glucose, for example [l|-Q. 

Nutrient consumption via facilitated diffusion is ubiq- 
uitous in nature and occurs in a variety of contexts, in- 
cluding oxygen consumption by human tissue cells, cal- 
cium uptake by intestinal cells, bacterial absorption of 
various sugars, etc 0-|g. In many of these examples, cells 
grow in clusters and shield each other from the available 
nutrients. In Fig. [TJ we see confocal microscope images 
of yeast colony cross sections in which yeast cells marked 
green arc growing due to an abundance of glucose. Nu- 
trient (glucose) shielding in larger colonies (Colonies 2 
and 3) prevents cell growth in the colony interiors (red 
regions in Fig. [T]) and only an outer shell of thickness £ 
is able to grow (illustrated for Colony 3 in Fig. [1]). 

Typical bacterial and yeast cell colonies (e.g., those in 
Fig. [T|) arc dense cell clusters with large cell packing frac- 
tions (j) > 0.5 [3, Q- In this paper, we characterize both 
densely and sparsely packed cell clusters. Dilute clusters 
with 4> < O.l can be constructed artificially, by suspend- 
ing microbial cells in a gelatinous matrix. Immobilized 
microbial cells have a very wide range of industrial and 
environmental applications (see Ref. [9| for a review). For 
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FIG. 1: (Color online) Confocal microscope images of cross- 
sections through the bottom of three budding yeast colonies 
(scale bars represent 100 fim). The red (darker) color is the 
constitutive expression of a protein, whose level is largely in- 
dependent of growth rate, in all cells and the green (lighter) 
color is ribosomal protein expression, indicating growth. 
Colony 1: 0.5 mM glucose, 43 h after inoculation; Colony 
2: 1.5 mM glucose, 47 h after inoculation; Colony 3: 4.5 niM 
glucose, 56 h after inoculation. Small colonies such as Colony 
1 receive enough nutrients for all cells to grow. Colonies 2 
and 3 are larger and growth occurs only in an outer shell of 
thickness I (illustrated for Colony 3). The red cells in the in- 
terior are shielded from the nutrients. We assume the colonies 
have a spherical cap shape (see Sec. lIVI for a discussion of the 
shape and experimental details). 



example, gelatinous beads are seeded with yeast cells and 
used in reactors to produce ethanol [l^l- To better un- 
derstand the growth dynamics and physical properties of 
these systems, it is important to characterize the nutrient 
transport properties of cell clusters as a function of both 
single cell nutrient uptake kinetics and the geometry of 
specific cell packings. 

A nutrient concentration ^p in some medium, such as 



water or gel, with a constant diffusion coefficient Dq 
obeys the diffusion equation 



dt^ - i^oVV- 



(1) 



In the steady state, the left hand side of Eq.[T] vanishes, 
and the equation reduces to Laplace's equation. This is 
an important and wcU-studicd equation in electrostatics, 
as it is the equation for the electrostatic potential $(r) = 
'0(r) in regions of space without any charges. 

In the case of nutrient diffusion, Eq.[T]must be satisfied 
everywhere in the medium outside of the cell. However, 
to completely solve Eq. [1] in the steady state, we must 
specify boundary conditions. One natural boundary con- 
dition is to set the concentration field at infinity to some 
constant value, i.e., ipilr] — > oo) — > ipoo, corresponding to 
a large nutrient bath with a uniform concentration ^oo- 
We also need boundary conditions on each cell surface. 
For example, if the cell is a perfect nutrient absorber, 
then for all points r on the cell surface S, the nutrient 
concentration vanishes, i.e. ■!/'(r) = 0. In the electrostatic 
analogy, this condition means zero electrostatic potential 
on every cell surface, i.e., each cell is a perfect grounded 
conductor. Conversely, if the cell does not absorb any 
nutrient, (i.e. it is a perfect refiector), then Fick's first 
law of diffusion tells us that the derivative of the con- 
centration along a direction n perpendicular to the cell 
surface must vanish. More precisely, the local nutrient 
flux density J(r) into the cell at some point r G S' on the 
surface satisfies 



J(r) := Don • VV'(r) 



(2) 



so that J(r) = for all r G 5 imphes n • V'0(r)|j.gg = 0. 
In the electrostatic analogy, this would correspond to a 
perfect insulator with no surface charge, with a vanishing 
normal electric field. Of course, living cells are neither 
perfect absorbers nor perfect reflectors. A more realistic 
boundary condition interpolates between these two ideal 
cases. 

A boundary condition on the cell can be derived from 
a more nficroscopic model of the nutrient transporters. 
For example, Berg and Purcell modeled transporters as 
small perfectly absorbing disks on the surface of an oth- 
erwise reflecting cell [ll|, [ij] ■ They showed that the cell 
requires very few transporters to act as an effectively per- 
fect absorber: A cell with as little as a 10"''' fraction of 
its surface covered by transporters takes in half the nutri- 
ent flux of a perfect absorber! Zwanzig and Szabo later 
extended this result to include the effects of transporter 
interactions and partially absorbing transporters |13l.ll4|. 
They showed that a homogeneous and partially absorb- 
ing cell surface model captures the average effect of all 
the transporters. As discussed below, in many cases of 
biological interest, the cell cannot be treated as a perfect 
absorber. The same partially absorbing boundary con- 
dition used by Zwanzig and Szabo will be derived in a 
different way in the next section. 

Although Eq.[l]is easily solved in the steady state for a 
single, spherical cell with the appropriate boundary con- 



ditions [ll|, [l^ , the complicated arrangement of cells in a 
typical mult i- cellular system, such as a yeast cell colony, 
implies a complex boundary condition that makes an ex- 
act solution intractable - one would have to constrain 
■(/'(r) and its normal derivative on a highly irregular ob- 
ject, like the surface of a cluster of grapes. In this paper, 
we explore an "effective medium" approximation to the 
exact solution of this problem. 

Effective medium theory treats a cluster of cells, or 
nutrient sinks, as a region with uniform effective nutri- 
ent transport properties (such as an effective diffusivity 
and nutrient absorption constant) that depend on the ar- 
rangement of cells in the cluster and the individual cell 
nutrient absorption properties. A key feature of the ef- 
fective medium theory is that these effective transport 
properties are derived in a self-consistent way. These the- 
ories have been used to calculate many effective proper- 
ties in heterogeneous systems such as conductivity, elas- 
ticity, and reaction rates (see [l5|, [lg| and Chapter 18 in 
[rn | for reviews). 

The paper is organized as follows: We develop our the- 
oretical model for nutrient uptake in single cells, dilute 
cell clusters, and dense clusters in Sec. |IT1 In Sec. IIIII we 
compare our analytic results with numerical solutions of 
Eq. [1] in the steady-state for clusters with hundreds of 
partially absorbing cells. We discuss experimental tests 
of our model in Sec. llVl and provide concluding remarks 
in Sec. IVl 



II. 



THEORETICAL MODEL 



We now discuss how to couple single cell nutrient up- 
take kinetics to the nutrient uptake behavior of an entire 
cluster or colony of cells via effective medium theory. In 
what follows we assume the cells are all identical and 
spherical. Although our experimental model is the bud- 
ding yeast cell, the theoretical treatment is quite general 
and can be adapted to any cell cluster that passively ab- 
sorbs nutrients that reach it via diffusion. 



A. Single Cell Nutrient Uptake 

We flrst review nutrient uptake by a single cell, with 
nutrient transporters operating via facilitated diffusion 
and following Michaelis-Menten kinetics. We smooth out 
the effect of a discrete set of transporters within the cell 
wall and study a radially symmetric model of nutrient 
uptake. A Michaelis-Menten model then means that the 
total nutrient current Icdi into a cell is related to the 
ambient nutrient concentration ipoo via 



-^ccll — 



V'o 



'■max v^oo 



Kjn + V'oo ' 



(3) 



where /max (sometimes called Knax in the literature) is 
the saturating nutrient current into the cell as i/'oo ~^ oo 



and K„i is the concentration at which In. 



t/2, as 





FIG. 2: (Color online) A plot of the nutrient current Jcdi 
(solid red line) into a cell as a function of the ambient nutrient 
concentration ?/joo, which follows Michaelis-Menten kinetics. 
The maximum current /max and the kinetics parameter Km 
are also shown. 



shown in Fig. [5] Even if a cell does not obey these kinet- 
ics for all ij^oo , one can often define a range of concentra- 
tions "000 characterized by effective kinetic parameters 
^max and Km, a characterization often used in experi- 
mental studies of nutrient uptake. Many studies infer an 
effective /max and Km from the measured amount of nu- 
trient consumed by a dilute suspension of cells 0, [iM \lM ■ 

The parameters /max and K„i are determined by the 
microscopic kinetics of each individual cell transporter 
and the density of these transporters on the cell surface. 
A third contribution arises from the structure of the cell 
wall, which influences the rate of transport of the nutrient 
into the cell. We can find /max and Km by appealing to 
a microscopic model of the transporters. Berg and Pur- 
cell, for example, modeled each transporter as a small, 
perfectly absorbing disk on the surface of a perfectly re- 
flecting cell [III, |T2] . Zwanzig showed that the large scale 
physics of this model are approximated by a partially ab- 
sorbing boundary condition on the surface of the cell [1^ ■ 
We will now derive this boundary condition in another 
way and connect /max and Km to our model. 

To model the cumulative effect of many microscopic 
details, such as transporter kinetics, the cell wall, etc., 
on the nutrient uptake we introduce a potential energy 
barrier U{r). Once a nutrient molecule surmounts this 
barrier, it gets absorbed by the cell. The nutrient con- 
centration ■ip{r) obeys the Fokker-Planck equation 



dt^ = DoV 



vv^- 



ksT 



iVU)'ip 



(4) 



where fcs is the Boltzmann constant and T is the tem- 
perature of the nutrient solution [2^ . Simple diffusion is 
recovered when the potential is constant. For simplicity, 
let's suppose that the nutrient must overcome a radially 
symmetric potential barrier U{r) that has a rectangular 
"lip" of height uq at r = a and with width w = a — a' 



FIG. 3: (Color online) A radially symmetric model of nutri- 
ent uptake in a single cell. We specify a steady concentration 
?/)oo far away from the cell. The cell is centered at the ori- 
gin and has radius a. The steady state concentration profile 
^Q{r) (black solid line) is calculated for a rectangular poten- 
tial barrier (red dashed line) of height no (designed to model 
the complex uptake dynamics of the cell wall) and a perfectly 
absorbing nutrient sink at r = a'. The concentration profile 
approaches tpao at large r and exhibits a jump discontinuity 
at r = a. The nutrient currents are constrained to be contin- 
uous. 



(see Fig. ig. Thus, 

U{r) = 



lio S{a — r) r > a' 
r <a\ 



(5) 



where d{x) is the step function {0{x) = 1, a; > and 
d{x) = otherwise) and uq is the height of the barrier. 
Wc assume perfect absorption at a' and a constant nu- 
trient concentration infinitely far away. 



^l;{r = a',t) =0 
ip{r — >■ oo, t) = ipa 



(6) 



To determine the nutrient flux into the cell, we solve 
Eq.|3]for the steady state profile ipQ{r). In addition to the 
boundary conditions (Eq. [5]), we ensure the continuity of 
the nutrient flux at r = a via the "jump conditions" at 
r = a, as discussed in Ref. [21|. The resulting concen- 
tration profile reads 



V'o 



aa'ipo 



1 



Mr) = 



l\gUo/ki 



a{r - a')-)pa 



(a' + (a- a')e^o/kBT'^ r 



r > a 



a' < r < a, 



(7) 



with tpo{r) = for r < a'. The shape of the solution 
4>o{i^) is shown in Fig. [31 

Let's now consider narrow potential barriers relative 
to the cell radius (w = |a — a'| <C a). Then, from Eq. [71 
we find the concentration gradient just outside the cell 
surface: 



SrV'O 



■ exp 



Up 

ksT 



Ma^) = Ki^o(a+), (8) 



where the + superscript indicates that we take the hmit 
r ^ a from outside the celL Eq. [5] reveals that the gra- 
dient of ■00 normal to the cell surface is proportional to 
'0o('') just outside. Notice that k — >■ oo when w — >■ (we 
also let uq — > 0), so that tpQ{a'^) — > at the cell surface to 
keep the flux finite. Thus, the cell is perfectly absorbing 
within our model if there is no potential barrier. Simi- 
larly, for a very large barrier [uq — > oo, ui finite), we have 
K — )• so that there is no flux of nutrient into the cell 
and dri^oiO''^) -^ 0, signifying a perfect reflector. 

The proportionality between a field and its gradient at 
a boundary is called a radiation boundary condition in 
the physics literature and can be derived quite generally 
[221 ■ This boundary condition is a natural coarse-grained 
description of the Berg and Purcell model of transporters 
as absorbing disks. Zwanzig and Szabo [ij, [ij] have 
used the radiation boundary condition to successfully 
model the physics of both perfectly and partially ab- 
sorbing disks on scales larger than the disk spacing, thus 
confirming our expectation that the coarse-grained nutri- 
ent uptake can be modeled by the ubiquitous radiation 
boundary condition with an appropriate choice of k. 

The absorptive strength of the cell can be parameter- 
ized by the dimensionless number v = na, where a is 
the cell radius. In chemical engineering, v is sometimes 
referred to as a Sherwood number [23|. If Pabs is the 
probability that a nutrient particle at the cell surface will 
be absorbed by the cell (instead of escaping to infinity), 
then first-passage techniques from probability theory [24 1 
lead to 



-'abs 



1 



(9) 



Thus, ly = Ka <S^ 1 indicates poor nutrient absorption 
while ly ^ 1 indicates a good absorber. Note that at i/ = 
1. the nutrient has equal probability of being absorbed 
at the cell surface or escaping to infinity. 

We now connect v with the measurable biological pa- 
rameters /max and Km- Recall that the nutrient flux into 
a single cell is related to the ambient nutrient concentra- 
tion -000 via the Michaelis-Menten relation Eq. [3] Sup- 
pose for now that the cells are well-separated, so the nu- 
trient uptake /coll of any given cell is independent of the 
others (i.e., there is no nutrient shielding), and that the 
ambient nutrient concentration ipo^ is held constant, so 
that the parameters /max and Km at each cell can assume 
their steady-state values. We assume the nutrient solu- 
tion experiences no macroscopic flows, such as convection 
currents, that would bias the isotropic absorption kinet- 
ics of the cell. The nutrient concentration ip{r = |r|) then 
satisfies V^V'(r) = in the steady state, with the bound- 
ary conditions ip{r — > cx)) = -000 and fi • V^/jj^ = Kf^ls at 
each ceU surface, as discussed above. 

Upon inserting Eqs. [7| and [Sj into Pick's first law [2J| 
(see also Eq. [5|) , we find the steady-state nutrient current 
into an individual cell, 

/cell = / /?on •V0(r)a dfi « — , (10) 

Js ^ + ^ 



where we integrate over the surface of the cell S (so that 
dil = sin 9 dd d(j) in spherical coordinates) and again as- 
sume w — a~ a' ^ a. Comparison of Eq. [TOj with Eq. [3| 
leads to 1/ as a function of biological parameters. In the 
limit of low ambient nutrient concentration (V'oo ^ Km), 
we have 

(11) 



AiraDoKm - /max 

It is also possible to define an effective v for a reflecting 
spherical cell uniformly covered by identical, partially ab- 
sorbing disks with radius fldisk and absorption parameter 
Kdisk- In this case, using the boundary condition Eq.[5|on 
each disk surface, Zwanzig and Szabo flnd that the effec- 
tive parameter v for the entire cell (for Odisk'^disk "^ 4/7r) 
is I' = A^diskadisk^disk/4a, where A'disk is the number of 
disks on the cell surface (see Ref. [l4|). 

A rough estimate of ly for glucose uptake by a 
S. cerevisiae cell follows from values for /max (4.2 x 
10"^ molecules/sec). Km (7.4 mM), a (2um), and Dq (670 
/im^/sec) found in the literature 0, l2a,l2y|- The /max is 
particulary difficult to estimate as the nutrient uptake 
rate in experiment is calculated per gram of dried yeast 
taken out of a liquid culture. To get an uptake rate per 
cell, we estimate that a yeast cell has a 2 pg dry weight 
[27| . We flnd that these yeast cells are in fact very poor 

absorbers with Pabs,yoast ~ I^ycast ~ 6 X 10""* -- 10"^ 

within an order of magnitude. Gram-negative bacterial 
cells are quite different! Again using literature values for 
/max (2 X 10^ molecules/sec). Km (1 A^M), and a (0.5 
/im) [iSJ . |28| for a single Escherichia coli cell, we find 

that Pabs,gram-n. ~ J^gram-n. ~ 0.09 -- 10"^ ThuS, this 

bacterium is ~ 10^ times more absorbent than the yeast 
cell: This striking difference has profound biological im- 
plications since, as we will show in the next sections, the 
parameter u greatly influences the growth of a cell colony. 

The large disparity in v values may be due to the 
thicker cell walls of S. cerevisiae, compared to gram- 
negative bacteria like E. coli. The presence of a cell wall 
can have two effects. First, the diffusion coefficient of 
the transported nutrients may be lower in the cell wall 
medium than it is in the bulk solution. Second, the ab- 
sorbing surface lies at the plasma membrane, not the 
surface of the wall. The second effect implies that the 
absorbing surface is at a distance a — w from the cell 
center (rather than a), where a is the cell radius and w 
the thickness of the wall. So, even if the diffusion coef- 
ficient inside the wall is the same as in the bulk and we 
have a perfectly absorbing plasma membrane, the effec- 
tive value of z^ is 1/ = a/w—\. This argument is consistent 
with the measured glucose uptake kinetics for the gram- 
positive bacterium Luconostoc mesenteroides, which has 
a cell wall thicker than E. coli and thinner than S. cere- 
visiae. We find Pabs,giamp. ~ t'gramp. « 5 X 10~^ from 
the literature values of the parameters [l^ [23| . Of course, 
other factors apart from the cell wall thickness could be 
relevant. 

The experimental results discussed here for single cells 
tell us that it is important to consider a large range of the 
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FIG. 4: (Color online) Schematic of the effective medium ap- 
proximation, which replaces A^ orange cells in a spherical clus- 
ter with radius b (enclosed by the dashed lines) in (a) by a 
homogeneous, attenuating medium in (6) shown in dark and 
light green. The light green rim of width £ illustrates the sec- 
tion of the cluster receiving enough nutrients to grow, xpao 
is the limiting nutrient concentration in the variably shaded 
red region, outside the cluster. To analyze this problem, we 
use the coordinate system in (c), with an origin at O, cell 
centers located at {R^}^!, and vectors Vi pointing from the 
cell center to the surface of the cell. 



parameter v = kg: v can range over at least two orders of 



magnitude (10 



< V 



< 10 ^) for yeast, gram-positive, 



and gram-negative bacteria. Thus, neither perfectly re- 
flecting nor perfectly absorbing boundary conditions are 
relevant for nutrient uptake in many cell populations. In- 
stead, we develop a theory for the nutrient absorption by 
a cluster of cells with arbitrary i/. This is the subject of 
the next section. 



B. Nutrient Uptake in Cell Clusters 

We now consider nutrient absorption by a cluster or 
colony of cells. In general, this is a very complicated 
problem involving solving the diffusion equation in the in- 
terstitial area of the cluster while making sure the bound- 
ary condition (Eq. [5]) is satisfied at each cell surface. Al- 
though analytical results are possible for a single cell, we 
must resort to numerical solutions and approximations 
when dealing with a cluster. A typical 1 mm diameter 
yeast cell colony contains over 10^ cells, making the exact 
solution for nutrient uptake by such a colony intractable 
even numerically. 

To model nutrient uptake (and nutrient shielding) in 
a group of cells, we again consider the steady-state dif- 
fusion equation with boundary conditions provided by a 
disorderly cluster of N identical spherical cells all with 
radius a and no overlaps (see Fig. 01). Let the cells be 



located inside a spherical region of radius b ^ a with 
centers at positions R°, where i = 1, . . . ,N. We also em- 
ploy local spherical coordinates at each cell to write vec- 
tors Ti = (a, f2i) pointing at positions on the cell surface, 
where Qi = {9i,(j)i) is a pair of polar angles that spec- 
ifics the direction of r^ relative to the center of the i*^ 
cell (see Fig. 111(c)). Consider a particular configuration 
of these N cells and denote by ^^(r) the exact solution 
to the steady-state concentration field in the interstitial 
region for this particular configuration. 

Upon assuming a time-independent steady-state, we 
modify Eq.[T]to include a set of Lagrange multiplier func- 
tions (7i{r) defined on the cell walls, 

N 
i—l 

We also set R^ = R° -I- Tj (see Fig. Hfc)) and choose r 
to be inside the cell cluster. The Lagrange multiplier 
functions {ai{ri)} will be chosen to satisfy the radia- 
tion boundary condition at each cell surface (the func- 
tions {o'i(ri)} would be charge densities in electrostatics). 
These boundary conditions (given Eq. [5] with v ^ na for 
a single cell) read, for all i and R^, 



*(Rz) 



r, • V*(R,) = Q,4'(R.), 



(13) 



where Qi = i^^^ r-i • V is a convenient gradient operator 
used in the detailed effective medium calculation in Ap- 
pendix [Xj The source function s(r) allows us to incorpo- 
rate additional boundary conditions on the concentration 
field. 

We now average over all possible cell configurations 
(consistent with excluded volume interactions between 
cells) to obtain the average transport properties of the 
nutrient in a cell cluster. Upon averaging both sides of 
Eq. [H we find 



/"s(r'-r)V(r')dr' + s(r), 



(14) 
(15) 



where ■0(r) = (^^(r)). The bracket average is an ensem- 
ble average over cell configurations and I](r) is a linear 
response function describes how the cells deform the con- 
centration field. 

The linear response approximation, justified here by 
comparisons with simulations, is only valid away from 
the cluster edges and for sufficiently small concentration 
field deformations. In addition, this approximation often 
breaks down for time-dependent diffusion [30| because 
the transient diffusive dynamics are dominated by slowly 
decaying modes due to large voids inside the cell clus- 



ter |3l| . Ref. [32l | uses a more microscopic description of 



nutrient diffusion to examine the validity of the linear re- 
sponse approximation in more detail. We will forgo these 



complications here and exploit the linear approximation 
above, checking our assumptions using experiments and 
simulations, as had been done for the physics of fluores- 
cence quenching [33|. 

An exact evaluation of S(r) involves an ensemble aver- 
age denoted by brackets in Eq. [TH This average requires 
the full probability distribution P({R^}^]^) of observing 
N cells with centers {R^jj^j^. Unfortunately, an exact 
solution obtained in this way would require knowledge of 
all of the correlations between the cell positions, which 
may not be experimentally accessible. We will assume for 
now that this distribution is known. Later, the effective 
medium theory developed in Sec. Ill El will approximate 
E(r) in a self-consistent way using just the one and two 
cell center distributions. 

What happens to the configurationally averaged so- 
lution ijj{r) of Eq. [15] over distances large compared to 
the size of a single cell? Specifically, how does the cell 
colony absorb nutrients on average, as if it were the ho- 
mogeneous medium illustrated in Fig. |31[6)? To answer 
this question, we perform a gradient expansion of Eq. 1151 
As discussed in Ref. [32|, such a gradient expansion ne- 
glects intrinsically non-local contributions to I](r) due 
to the fluctuations in the concentration field that can- 
not be averaged over large distances. The mean field 
approximation used here then uses the resulting trans- 
port coefficients to describe the absorptive properties of 
the cluster. Simulations have shown that this approach 
correctly models the physics in related systems |34| . 

Eq. [TS] in Fourier space reads 



- Doq^^j{ci) = S(q)VXq) + s(q). 



(16) 



where we have applied the convolution theorem and all 
the functions arc now their Fourier transformed functions 
of the 3d wave-vector q. On average, the cells in the 
cluster should be distributed isotropically and S(q) can 
only depend on q = |q|. Expanding E around q = gives 
us the desired gradient expansion: 



-I?ogV(q) = S((z - O)^(q) 



S"(q = 0)g2V-(q) 



+ 0{q^) + s{q) 
= kiPiq) + 5Dq^^iq)+Oiq^) + siq), (17) 

where we have identified an absorptive term kip and a cor- 
rection to the diffusion term 6Dq^t/j. We neglect higher 
powers of g (i.e., higher order derivatives of the concen- 
tration field) in our coarse-grained reaction-diffusion de- 
scription of nutrient transport in a large colony of cells. 
Upon returning to real space, we find the desired macro- 
scopic transport equation for the configurationally aver- 
aged nutrient concentration: 



DV^iP{r) - fc^(r) = s(r) 



(18) 



where D = Dq + 5D. The coefficients D and k charac- 
terize the macroscopic diffusion and absorption, respec- 
tively, of the nutrient in the cell colony. The crucial 



question is how the effective diffusion constant D and 
absorption k depend on parameters such as v and the 
cell volume fraction (j). 

In general, wc expect that k will be a positive, increas- 
ing function of the cell volume fraction 0, since the cluster 
becomes more absorbent as wc introduce more cells. It 
will also increase with v = Ka as each cell will absorb 
more nutrients as v increases. The sign of 5D is more 
subtle, because of two competing factors: The confine- 
ment of the nutrient by the cells in the interstitial space 
will decrease the effective diffusion constant D, while the 
nutrient gradients in random directions induced by nu- 
trient uptake by cells will increase D. Since v controls 
the amount of absorption, we expect that (for small 0) 
SD < for I' '^ 1 and 5L» > for i/ > 1. This is 
confirmed by our effective medium calculation of D in 
SecHHl 

Wc assumed that the radii of the cells in the cluster 
arc monodisperse. This approximation might not be re- 
alistic in cell clusters at different stages of their cell cycle. 
However, in both polydisperse and monodisperse cases, 
we can still define a cell volume fraction </>. The effective 
transport coefficients D and k will depend on this vol- 
ume fraction and the polydispersivity of the radii. For 
a fixed 0, polydisperse cell radii will decrease the total 
cell surface area and, consequently, decrease the nutrient 
absorption k. For example, the ratio of the total cell sur- 
face area S'poiy of a cluster of cells with a Gaussian radius 
distribution (with average (a) and variance aa <C (a)) to 
the total surface area S'mono for a cluster of cells with the 
same radius (a) is given by [17[ 
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(19) 



Thus, a theory of monodisperse cells is reasonable pro- 
vided (a) > VSaa- Also, we don't expect the cor- 
rection to vary much with u since the main effect of 
polydispersivity seems to be from the surface area de- 
crease. Finally, computer simulations for two discrete 
sphere sizes reveal that the ratio fcbidisp/^mono of absorp- 
tion coefficients for (perfectly absorbing) cells satisfies 

"^bidisp/fcmono ~ (>-'bidisp/'-'mono) [351. loDJj OVCr a Wide 

range of packing fractions. Thus, it may be possible 
to approximate polydispersity in our theory by reducing 
the absorption coefficient k by such a geometric factor. 
However, unless stated otherwise, we henceforth ignore 
this complication and instead consider cell clusters with 
monodisperse radii. 



C. The Macroscopic Screening Length ^ 

An important nutrient screening length associated 
with Eq. [THl for cell clusters is 



(20) 



We now use effective medium theory to calculate ^ and 
relate this length to cell configurations and single cell 
nutrient uptake kinetics. 

Consider first a spherical cell cluster (or "colony") of 
radius b in which ip satisfies the effective medium re- 
sult Eq. [m Assume as well that V^V = outside of 
this spherical region. Then, with the boundary condi- 
tion tplr — > oo) = -000 and continuity of the nutrient 
current at r = 5, wc find 
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(21) 

where T{x) = tanh(a;)/a;. When b ^ ^, the concentra- 
tion field 0(r) Ri °^^^ exp[(r — b)/^] near the colony 
surface, and thus decays exponentially as we move fur- 
ther into the interior. Thus, ^ is a nutrient penetration 
depth. 

Eq. ED leads to the total flux /duster into the cell cluster 
in the same way as the single cell discussed above. We 
will see in the next two sections that a good approxima- 
tion is D Ki Dq. Then, the nutrient flux into the cluster 



/cluster = 47r_DoV'o 



& — ^ tanh 



(22) 



Thus, when ^ ^ 6 so the nutrient does not penetrate 
far into the cell colony, /duster ~ 47r/)o'0oo& and the en- 
tire colony acts as if it is a single, perfectly absorbing 
sphere with radius b. Conversely, if ^ ^ 6, then nutrient 
penetrates deep into the colony. In this limit we have 
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(23) 



so that the total nutrient absorption is down relative to 
a perfectly absorbing colony by a factor oi S^/3 ^ 1. 
/duster now scales with colony volume (oc b^), since each 
cell in the colony contributes to the nutrient uptake. 

A related biologically relevant parameter is the thick- 
ness £ of the outer shell of actively growing cells, as il- 
lustrated in Fig. [Hand by the green band in Fig. \M,b): 
Suppose that cells require some minimum concentration 
V'min of nutrient in order to grow (or exhibit some level of 
growth-coupled fluorescence) . We can estimate £ by flnd- 
ing the location r^m^ inside a spherical colony such that 
'^(''min) = "^min (sec Eq. [2T|) and setting £ = b-r„,in- For 
b ^ ^ and D « Dq, this length is related to b, ^, ipn-iin 
and ■i/'oo via (see Fig. [S]) 
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where W is the Lambert- Vl^ function [37|. The thickness 
£ can also be determined experimentally using a fluores- 
cence reporter, as shown in Fig. [TJ Experimental esti- 
mates of £ in yeast cell colonies are discussed in Sec. IIVI 
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FIG. 5: (Color online) Plot of tfie thickness £ of actively grow- 
ing cells in a spherical cluster of size & as a function of the 
nutrient bath concentration ^oo for various values of the pen- 
etration depth ^. 



D. The Dilute Cell Limit 

As mentioned in SeclH dilute clusters of microbial cells 
can be realized in experiment by embedding the cells in 
a matrix. The initial cell density in such a matrix can be 
tuned over a wide range and is often quite small. More- 
over, since the matrix usually has a negligible influence on 
the diffusion of small molecules like glucose [S^ , we now 
develop an effective medium theory for clusters with cell 
densities low enough to neglect interactions. This theory 
will serve as an important and instructive limiting case. 

We approximate the nutrient flux into each cell with 
the single cell result Eq. [TUl replacing tpoo with the lo- 
cal value ■0(1") I where r is the location of the cell. This 
approximation does not take into account the flnite size 
of our cells or the deformation of the concentration field 
^(r) around each cell. Thus, we cannot use this dilute 
limit to find a correction to the diffusion constant D. 
However, as we will see in the next section, we expect 
this to be a very small correction, especially for small 
packing fractions (j). 

To approximate the nutrient absorption k, we assume 
that each cell independently contributes to a nutrient flux 
per unit volume 
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where n is the cell number density which we relate to the 
cell volume fraction </> = | na^n. We conclude that the 
macroscopic transport coefficients in this dilute approxi- 
mation are 
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More sophisticated theories can find a correction to D 
that is linear in (j) and higher order corrections to k (see 




FIG. 6: (Color online) Plot of the screening length ^d in the 
small limit as a function of v for various values of </) on a 
log-log scale. ^ crosses over to the perfectly absorbing limit 
at !/ ~ 1. The arrows Y and B denote, respectively, typical 
values of v for baker's yeast and bacteria. 



Rcf . [lg| for a review) . Eq. [21] implies that the screening 
length ^ = ^ D jk in the dilute limit is given by 
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Fig.|6]shows ^d as a function of the absorption parameter 
i> = Ka: the screening length is a strong function of v for 
v < 1 and crosses over to a perfectly absorbing regime 
when 1/ > 1 . If the cells in a dilute cluster at low volume 
fraction have polydisperse radii in a Gaussian distribu- 
tion of cell radii with average (a) and standard deviation 
fa <C (a) , Eq. [27] becomes 
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As expected, the reduced total cell surface area in the 
polydisperse cluster leads to a larger screening length. 
Note that the ratio S,d, poiy/Cd is independent of v and 
in this limit. 

For yeast cehs with « 10"^ and u « 10"^, Eq. [27] 
predicts ^ ~ 60a. Thus, glucose penetrates far into the 
yeast cell cluster when </) w 0.1 (e.g., for yeast embedded 
in gel) , allowing a substantial fraction of the yeast pop- 
ulation to grow. Conversely, we expect that ^ sa 6a for a 
bacterial colony [v k, 10~^) at w 0.1 volume fraction. 
We clearly need to go beyond the dilute limit and de- 
termine the dependence of ^ on and v more carefully: 
Yeast and other cell colonies rarely grow at low volume 
fractions; cells typically clump together and pack them- 
selves in an amorphous structure with a volume fraction 
approaching that of the random close packing density, 

(^w 0.6 -0.7 m. 



E. Effective Medium Theory for Dense Cell 
Clusters 



For (/) > 0.1, the dilute approximation breaks down, 
and we must solve Eq. [T2] for \l'(r) more exactly. 
Upon defining the diffusive Green's function Go(r) ~ 
(47r£'o'')~^, we can rewrite Eq. [T2] as an integral equa- 
tion, 

*(r)= /'dr'Go(r-r') 

- N 

X ^ /daa,(r,)(5(r'-R0 + s(r') , (29) 



where we integrate over all positions r' (dr = d^r). In 
principle, we could choose the functions ai{ri) to enforce 
the radiation boundary condition at each cell surface 
(Eq. [13]), and then average over the cell positions to ob- 
tain S(r) in Eq.[T5] However, as discussed in Section lllBl 
we do not have access to the full probability distribution 
P({R"},j^]^) necessary to do the averaging. An alter- 
native approach exploits an approximate solution which 
accounts for each cell independently to first order, the 
effects of pairs of cells at second order, etc. This expan- 
sion requires knowledge of the probability distributions 
of a single cell position, a pair of cells, a triplet, etc. 
Such scattering expansions are plagued with divergences 
due to the long-range nature of the diffusive interaction 
between the cells and must be treated via careful resum- 
mations [32, Eo, EH . Particular care is required, because 
the Green's function, or propagator, in this expansion is 
G'o(q) = (^09^)"^ in Fourier space, which is singular as 
g-^0. 

To bypass such complications, we imagine cells imbed- 
ded in an "effective medium" with effective transport 
properties, to be determined sclf-consistently. These 
transport properties include screening which cuts off the 
long range behavior of Go(r), and renders G'(q) finite as 
q — >■ 0. In particular, we assume a modified "effective 
medium" propagator given by 



Gs(q) 



i5og' + E(q)' 



(30) 



which describes the transport properties of the nutri- 
ent as it wanders through the homogeneous attenuating 
medium shown in Fig.[4l[6). A self-consistency condition 
on S(q) leads to a renormalized (and better behaved) 
scattering expansion. 

In Appendix \K[ we follow Cukier and Freed's effective 
medium calculation of S(q) [42] . However, our cells will 
be partially absorbing, which generalizes the perfectly 
absorbing case Cukier and Freed considered. The calcu- 
lation assumes that the cells have uniformly distributed 
centers in the cluster. Excluded volume interactions are 
included in an approximate way, by assuming that the 
centers R^ ■ of any pair of cells i and j are distributed 
according to the low density hard sphere pair distribution 



OS 

f 0.6 

u 

0.4 
0.2 


• • 


*" *fc .^ 


^— 0.001 
- • • Cl.l 

1 

•• 5 

•■ ion 


""'"■■•'■'■'.,. 






, , , 



0.0 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



FIG. 7: (Color online) Plot of the screening length ^ calcu- 
lated with the effective medium theory divided by the dilute 
limit result ^d as a function of the volume fraction <f) for var- 
ious values of v. The low density result becomes increasingly 
inaccurate at larger packing fractions. 



1.3 
1.2 
1.1 

D 1.0. 

Dq 0.9 

0.8 
0.7 
0,1 



.''l--''" ' 




i --- 0.1 

.... 1 

: 6 

100 


~* --i 



).0 



0.1 



0.2 



0.3 





0.4 



0.5 



0.6 



FIG. 8: (Color online) Plot of the diffusion coefficient D in 
the cell cluster calculated with the effective medium theory 
divided by the bare diffusion coefficient Dq in the absence of 
the cluster as a function of the volume fraction for various 
values of v. 
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where 9{x) is the step function and V is the cluster vol- 
ume. This generalization was also considered by Cukier 
[43| , but we find an important second order correction 
to his results from the pair distribution function Eq. [31] 
The details are contained in Appendix 1X1 

Wc find that the nutrient screening length in units of 
the cell radius a = ^/a satisfies 
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where ie{x) and ki{x) are the modified spherical Bessel 
functions of the first and second kind, respectively. For 
a specific v ~ Ka, Eq. [32]is a self-consistent equation for 
a = £,/a. It is soluble numerically (we used Newton's 
method of successive approximation) for ^ as a function 
of u and (f>. 

The ratio ^/^d, where ^d is the dilute limit formula 
Eq. [571 appears as a function of (j) for various v in Fig. [71 
Note that the effective medium theory always predicts 
^ < ^d, a plausible result since this more sophisticated 
theory allows for nutrient absorption mediated by re- 
peated interactions with individual cells. 

For the values of u typical of yeast clusters and gram 
negative bacteria {y ^ lO^'^ — 10"^), ^ is no more than 
20% different from S,d even for large (/) ^ 0.5. These 
results suggest that effective medium theory is a rea- 
sonable approximation for dense clusters of cells with 
V = na <^ 1, m. contrast to the more problematic per- 
fectly absorbing case i^ — )■ oo |60 |. 

The second order term in the Taylor series expansion 
of E(q) around q = yields an implicit equation for 



the effective diffusion coefficient D of the nutrient inside 
the cell cluster. The somewhat cumbersome equation 
is Eq. IA20I in Appendix \^ The solutions for various v 
are plotted in Fig. [8l as a function of 0. In all cases of 
biological interest, the corrections to the bare diffusion 
coefficient Dq in the absence of the cluster are relatively 
small, i.e., approximately 10-20% and never exceeding 
30%. 

In the biologically relevant regime v < 0.1. the cor- 
rection to the diffusion coefficient is close to the i/ = 
(perfectly reffecting) limiting case given by 
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(33) 



which is derived from Eq. IA201 As expected, the dif- 
fusion coefficient D is smaller than Dq when the cells 
are perfectly reflecting. Also, the structure of Eq. [331 is 
consistent with the modified Maxwell- Garnett theory for 
diffusion of particles in a medium with perfectly reffecting 
inclusions derived in Ref . 4J] . Eq. [331 is also consistent 



with experimental results for the diffusion coefficient of 
water in a suspension of spherical colloidal particles [43 ■ 
At high volume fractions, we expect the effective medium 
theory to start to break down and higher order terms in 
E(q) to contribute to transport within the colony. Thus, 
the curves in Fig. [H should be treated as approximations 
and specific features, such as the crossover between the 
1^ = and v = 0.1 lines, might not be physical. 

Wc now check the effective medium result against sim- 
ulations and experiments. 



III. NUMERICAL SIMULATIONS 

To obtain further insight into nutrient absorption, we 
now go beyond effective medium theory and solve the 
diffusion equation numerically in a densely randomly 
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packed cluster of spherical cells. We employ the fi- 
nite element solver discussed below to numerically solve 
the steady-state diffusion equation with the appropriate 
boundary conditions. Given the complicated nature of 
the problem, we can only test effective medium theory for 
clusters of up to 400 cells in this way. However, since the 
effective medium approximation is derived in the limit 
of an infinite number of cells, if it works for simulations 
with 100-400 cells, it should be even more reliable when 
we have even more cells, as in a typical growing yeast 
colony. 



A. Cell Cluster Simulation 

To check the analytic results, we solved the steady 
state diffusion equation exactly for cell clusters with hun- 
dreds of cells. The numerical solution was found with the 
COMSOL 3.5a finite element solver g^. A MATLAB 
program was written to input in the locations and radii 
of all the cells in the cluster. The COMSOL program 
included a computer-assisted design (CAD) feature that 
was then able to parse the MATLAB output and cre- 
ate a particular arrangement of spherical cells that de- 
fined our domain of interest. The coordinate list for the 
sphere cluster was created via a Bennett model, originally 
designed to quickly generate amorphous, dense random 
packings of identical spheres [43]. These arrangements 
approximate the disordered packing of cells observed in 
yeast colonies in the experiments. 

To prevent the cells from touching and creating singu- 
larities in the finite element mesh, identical spherical cells 
with radius a are placed at the sphere centers of a clus- 
ter generated by the Bennett model using spheres with 
larger radius a > a. This guarantees a gap of at least 
2 (a — a) between adjacent spheres. A high volume frac- 
tion (f) « 0.63 (corresponding to random close packing) 
is generated with a = 1 and a = 0.999. Smaller values 
of (j) are generated by decreasing a for a fixed a = 1 . For 
example, a = 0.9 results in a volume fraction (j) w 0.48 
(e.g., the cluster in Fig.[S|). 

A large bounding sphere concentric with the center of 
mass of the cell cluster allowed us to impose a constant 
nutrient concentration "at infinity": ^(|r| — >■ oo) = V'oo, 
to approximate a suspended cell cluster in an infinite 
nutrient bath buffered at concentration ipoo. We used 
COMSOL's "infinite element" option to efficiently sim- 
ulate by placing extremely large finite elements between 
the bounding sphere and the cell cluster. With bounding 
spheres of radii typically five times the cell cluster radius, 
this method yielded good results: Our calculated solution 
was unchanged when varying the bounding sphere radius 
from 5 to 10 times the cluster radius. 

Finally, we specified the dimensionless nutrient uptake 
parameter v = na for cells with identical radius a and ap- 
plied the boundary conditions given by Eq.[T3]on each cell 
surface. A particular cell cluster and the corresponding 
steady state nutrient concentration ^(r) near each cell 
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FIG. 9: (Color online) A simulated cluster of cells with u — 1 
and tpoo = 1. The color indicates the local nutrient concen- 
tration 'l'(r) near the cell surfaces. In this case, the 200 cells 
with radius a were confined to a cluster of radius b ^ 7.45a, 
yielding a volume fraction of (^ « 0.48. The radius of the 
large bounding sphere on which the nutrient concentration 
was fixed at tpoo was about 56a ~ 7.5b. To increase compu- 
tation efficiency, we used COMSOL's "infinite element" op- 
tion with spherical symmetry between radii 28a and 56a (see 

Sec.lmll. 



surface is shown in Fig. [HI In this case, 1/ = I, the packing 
fraction cj) w 0.48, and the cluster radius b sa 7.45a, with 
f/joo = 1. The bounding sphere had a radius of 56a and we 
inserted infinite elements at distances between 28a and 
56a. With these parameters, the effective medium theory 
predicts a screening length ^ w 0.72a. This is consistent 
with what we observe in Fig. [5J as the concentration de- 
cays by a factor of e over a distance coinparable to the 
cell radius a. 



B. Comparison with Theoretical Results 

To compare theoretical and simulation results, we 
tracked the decay of the radially averaged nutrient con- 
centration 'l>(r) = (47r)^^ J *(r) dfl into the center of the 
colony, and compared with the effective medium predic- 
tion ip{r) of Eq.[lTJ Fig. [TU] shows a semilogarithmic plot 
of r'iJj{r)/a'ipoo versus 6 = (b — r)/a, varying v = na over 
3 orders of magnitude. Eq. [21] predicts that this quan- 
tity decreases exponentially near the cell surface when 
1 <C b/a. Fig. [TU] shows how the nutrient concentration 
decreases as we move into the cluster. 

Fig. [TOl shows that, even for the large volume fraction 
(j) ~ 0.63, effective medium theory provides an excellent 
description, especially for small values of v. Thus, this 
theory is appropriate for modelling yeast colonies, which 
have a very small v ^ 10~^ value. The effective medium 
agreement should improve if computer resources allow 
more cells in the simulation, since it was designed to han- 
dle the limit where the cell number iV — >■ 00 and cluster 
volume y — >■ 00, with N/V fixed. Fig. [TT] shows the low 
V regime (poor nutrient absorbtion) in more detail. Al- 
though the differences between the simulation and effec- 
tive medium theory are now more evident, the absolute 
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FIG. 10; (Color online) The radially averaged concentra- 
tion r^(r)/a?/'oo (symbols) and the effective medium result 
rip{r)/atpoo (lines) for a cluster of 353 cells of radius a (with 
cluster radius b ~ 8.15a) as a function of J = (6 — r)/a on a 
log-linear plot. The concentrations are rescaled to highlight 
the exponential decay of the concentration into the bulk of 
the cluster. Cells occupy the region < 5 < 8.15 with the 
center at 5 « 8.15, as indicated by the dashed lines. The cell 
packing fraction is (^ f» 0.63. The bounding sphere has a ra- 
dius of 50a, with variable mesh elements inserted at distances 
between 25a and 50a. 
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FIG. 12: (Color online) Exponential decay of the rescaled 
nutrient concentration in the effective medium theory (lines) 
and in a simulation (symbols) for a more dilute 85 cell cluster 
with b — 7a and various v. The cell packing fraction is </!> « 
0.24. The cluster, shown above the graph, now occupies the 
region between the dotted lines, < 5 < 7. 
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case is obtained for the highly absorbing limit v > I, 
as well. We conclude that effective medium theory pro- 
vides a good description of spherically averaged nutrient 
uptake, for both weakly and highly absorbing cells, for 
volume fractions 6 < 0.6. 



IV. EXPERIMENTS 



FIG. 11; (Color online) The radially averaged concentra- 
tion '^{r)/tpao (symbols) and the effective medium solution 
'4>{r)/'4'oo (lines) as a function oi 5 — {b — r)/a for the same 
353 cell cluster described in Fig. 1101 The dotted lines indicate 
the cluster edge {S — 0) and center {S ~ 8.15). 



difference between the simulation and the theory for the 
concentration remains small. 

We also studied a more dilute 85 cell cluster, to check 
that effective medium theory is indeed accurate for lower 
volume fractions and an alternative cell configuration. 
Instead of using the Bennett model, we placed cells with 
centers in three concentric shells with radii 2a, 4a and 
6a. The cells within each shell were placed randomly, 
but their positions were adjusted to prevent the cells 
from touching and disrupting the finite element mesh. 
Fig. [12] reveals even better agreement between the simu- 
lation and effective medium theory for the lower volume 
fraction of « 0.24. Note that good agreement in this 



We now compare the effective medium theory pre- 
sented in Sec. |ll] to experimental results for glucose up- 
take in yeast cell colonies. Yeast cell aggregates are a par- 
ticularly interesting biological application because they 
form naturally in the wild and are a possible model for 
the initial emergence of multicellularity [i^. Nutrient 
uptake in single versus clumped cells is a crucial factor in 
this model system. Thus, accurately modelling nutrient 
transport properties in cell clusters is key to understand- 
ing their biological function and role in evolution. 

Since yeast cells are poor glucose absorbers, the char- 
acteristic nutrient penetration depth in a yeast colony 
is much larger than a cell radius and is easily visible in 
colony cross sections, as shown in Fig. [TJ Their poor 
absorption properties are also well-suited for effective 
medium modelling, since we expect that nutrient trans- 
port for small absorption parameter i/ <^ 1 is dominated 
by much longer length scales than the local colony geom- 
etry at the single cell scale. 
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(a) Confocal microscopy setup 




[b] Fluorescence of irdividual cells 
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FIG. 13: (Color online) (a) Confocal microscopy setup. A sin- 
gle yeast cell was inoculated into a well of a glass-bottomed 
96-well plate containing 200 /iL of yeast media and allowed 
to grow into a colony. The bottom of the colony was imaged 
using a spinning-disc confocal inverted microscope. (6) Cells 
were grown in shaken media, pipetted into an empty well, 
and imaged immediately (before they could divide and form 
colonies) with the setup shown in (a). At least 250 cells at 
each concentration were imaged, and the error bars are the 
95% confidence interval using the Student's t-test. Relative 
fluorescence intensity in (6) is deflned as ribosomal protein ex- 
pression divided by constitutive expression. Cells in 0.167 mM 
glucose continue to divide, but at less than half the growth 
rate of cells in 0.5 mM glucose. 



A. Experimental Setup 

We used confocal microscopy to measure the distribu- 
tion of growth in a yeast colony. To start the experiment, 
we filled each well of a 96-well, glass-bottomed microliter 
dish (Greiner Bio-One 655892) with 200 fiL of minimal, 
synthetic yeast media (recipe in Ref. ji^). Figure \W[a) 
shows the confocal setup and Fig.|T]shows representative 
images. We engineered a yeast strain (yJHK192) with 
two fluorescence protein (FP) reporters: (1) Cerulean, a 
stable CFP (cyan) variant, was constitutively expressed 
from the ACTl promoter; and (2) mCitrine, a YFP (yel- 
low) variant, was expressed from the ribosomal protein 
RPL7A promoter. The expression of RPL7A has been 
shown to correlate with growth rate [4^; the correla- 
tion between growth rate and glucose concentration has 
been well studied in yeast [50], and we verified that the 
YFP fluorescence correlated with glucose concentration 
in cells grown in shaken media (Figure IT^f fo)). mCitrine 
was destabilized by adding a CLN2 PEST-sequence to 
the C terminus [51j : a destabilized fluorescent protein 



was essential to see decreased expression of RPL7A in 
cells that had stopped growing. 

When grown on agar, small yeast cell colonies (less 
than 1 mm across) form a spherical cap with a contact an- 
gle 9 between the edge and the agar surface. The contact 
angle 9 increases with agar concentration and 9 k, 40° for 
2% agar [52l |. Unfortunately, in our experiment, the side 
view of the colonies was optically inaccessible and a pre- 
cise characterization of the colony morphology was not 
possible. However, since our colonics are very small, we 
expect that surface energy minimization is most impor- 
tant in determining the shape so that the colonies are 
approximately spherical caps and close to hemispherical. 
Also, we know the edges of the colonies were at least 
three cells thick (the maximum depth we were able to 
see by microscope). 

We grew colonies in three glucose concentrations: 0.5 
mM, 1.5 mM, and 4.5 mM. We inoculated one cell 
per well using a fluorescent activated cell sorter (MoFlo 
FACS, Beckman Coulter, Inc.); the cells had been grown 
to saturation in 1 mM glucose synthetic media. Inocu- 
lation was verified under a microscope, and those wells 
with a cell closest to the middle of the well (at least 3 
wells per glucose concentration) were selected for colony 
tracking. The plate was incubated without shaking at 
25 °C for 3 days. Each colony was imaged a day after 
inoculation and an additional 2-3 times over the next 2 
days. The radius of the largest imaged colony in each 
glucose concentration was less than half of the average 
colony size after two weeks of growth (610 m in 0.5 mM 
glucose, 920 m in 1.5 mM glucose, and 1170 m in 4.5 
mM glucose), indicating that the carbon source was not 
yet depleted. Images were taken with a 20X objective on 
a Nikon inverted Ti microscope with a Yokagawa spin- 
ning disc unit and an EM-CCD camera (Hamamatsu Im- 
agEM); CFP was excited with a 447 nm laser and YFP 
was excited with a 515 nm laser; exposure times for all 
images (including single-cell images) were 200 ms (CFP) 
and 1000 ms (YFP). All images were focused on the bot- 
tom layer of cells in the colony, and multiple, overlapping 
images were taken of colonies that exceeded the field of 
view. Three independent experiments were performed. 
Figures [T4| and [T5| show the combined results of all ex- 
periments. 



B. Data Analysis and Results 

Images were processed using the Fiji distribution of 
ImageJ 1531 . Images were converted to 8-bit, stitched to- 
gether [5J], and merged into a single RGB image. Con- 
trast was not adjusted during processing. Fluorescence 
intensity as a function of colony radius was measured 
using a custom script written in Python. The basic algo- 
rithm is as follows: (1) The constitutive (CFP) image was 
thresholded using Li's Minimum Cross Entropy thresh- 
olding method [55|. (2) Noise and cells not attached to 
the colony were removed by eroding and then dilating the 



13 



2.5 -f 

2.0 

1.5 

1.0 

0.5 

0.0 
2.5 

2.0 

1.5 

1.0 

0.5 

0.0 
2.5 

2-0 

1.5 

1.0 

0.5 

0.0 




Q.S 
mM 



■250 




1.5 
mM 




HwtaiiiMMMHMHM 



T 



T 



mM 



di 



1 — 

100 200 300 400 

stance from colony edge d = b r (/inn) 



FIG. 14: (Color online) Relative fluorescence intensity (ribo- 
somal protein expression divided by constitutive expression) 
as a function of distance d = h — r from the edge of the colony. 
Here b is the colony radius and r is the observation position. 
Each point is the average fluorescence of an 8 /xm thick ring 
whose edge is located an equal distance from the edge of the 
colony. The last point on each connected set of lines repre- 
sents the radius of the colony and larger radii reflect longer 
periods of incubation. Each line is one colony. Colonies were 
grown in 0.5 mM glucose (top), 1.5 mM glucose (middle), and 
4.5 mM glucose (bottom). Profiles from the three colonies 
shown in Fig. [T] are shown as orange triangles. 



binary image a total of three times with a four pixel di- 
ameter circular structuring element. (3) A series of mask 
rings was made by a series of morphological erosions and 
subtractions from the original image. The distance be- 
tween rings was approximately 8 /xm, or 2 cell diameters. 
(4) The masks were used to generate images of concentric 
rings using both the constitutive (CFP) and the growth- 
dependent expression (YFP) image. (5) The average flu- 
orescence of each ring was measured; the reported, rel- 
ative fluorescence is the average expression (YFP) fluo- 
rescence divided by the average constitutive (CFP) flu- 
orescence. Figure [14] shows the intensity as a function 
of distance from the edge of the colony; each point in 
the flgure is a measurement of fluorescence at one of the 
concentric rings, and a line joins the values from a single 
colony. 

In order to compare the experimental results to the 
effective medium theory, we assume that each colony 
is a hemisphere with a radius h equal to the radius of 
the bottom cross-section. Since the colonies are small 




100 200 300 400 

colony radius h (^m) 

FIG. 15: (Color online) Penetration depth ^ as a function 
of colony radius b. Penetration depth in the experiment 
(symbols) is the distance from the edge of the colony at 
which relative fluorescence drops below 0.4. The theoreti- 
cal results (solid lines) are calculated by numerically solving 
il>{r — b — I) = Vmin (Eq. I21|) for £, where we assume that 
^min ~ 0.25 mM local glucose concentration corresponds to a 
0.4 relative fluorescence level. As discussed in the text, this 
level is slightly higher than the level at which the cells stop 
growing and expressing the ribosomal protein. We also use a 
cell radius of a = 2 ^m, a cell packing fraction </> = 0.56, and 
!/ = 6x 10"". 



compared to the glass well dimensions and since the glu- 
cose molecules do not stick to the glass, we also assume 
that the hemispheres sit on an inflnite, reflecting sur- 
face. Hence, by reflection symmetry, the concentration 
proflle ■(/'(r) is given by the solution for spherical colonies 
(Eq. I^Tj) . where r is the distance from the colony cen- 
ter (measured above the inflnite plane). Of course, as 
discussed in the previous section, the colonies in the ex- 
periment could be spherical caps with a smaller contact 
angle 6 < 90° at the glass well bottom. A smaller con- 
tact angle corresponds to a more shallow colony with an 
increased nutrient penetration depth, due to nutrient dif- 
fusion from the top of the colony. However, we expect 
that the corrections due to smaller 9 to be small when 
the penetration depth is small compared to the colony 
radius. 

Some predictions of the effective medium theory are 
qualitatively conflrmed by the experimental data. The 
effective medium theory described in Sec. [H] predicts that 
for these dense cell colonies, the characteristic shielding 
length is 50 /im, or about 12 cell diameters. This is con- 
sistent with the fluorescence curves in Fig. [TD The the- 
ory also predicts that the glucose level at the outer edge 
of the colony decreases like 1/6, where b is the colony 
radius (see Eq. [?T|) . We see this effect in Fig. [T31 where 
shorter curves (corresponding to smaller colonies) have a 
higher fluorescence level at d = 0. However, the fluores- 
cence level has a complicated, non-linear relationship to 
the glucose concentration in the bulk medium shown (see 
Fig.HU^)) and the specific 1/6 scaling cannot be tested. 

Figure [TSl shows the penetration distance ^ as a func- 
tion of colony size; penetration distance is defined here as 
the distance from the edge of the colony at which the rel- 
ative fluorescence intensity falls to less than 0.4. This in- 



tensity is chosen to be well above the background fluores- 
cence level. For the smaller colonies, the fluorescence is 
higher than 0.4 throughout the colony. For these colonies, 
we say that the penetration depth is equal to the colony 
radius: £ = b. Eq. [21] is used to find a theoretical result 
for the penetration depth £ inside the colony. We assume 
that the relative fluorescence intensity of 0.4 corresponds 
to a local glucose concentration of V'min ~ 0.25 mM. 
Again, note that this level is slightly higher than the 
minimum glucose level required for growth. We estimate 
the packing fraction of cells inside a colony, cj) « 0.56, by 
looking at cell colonies growing on agar (data not shown). 
Then, assuming a cell radius a = 2 /im and absorption 
parameter ly — 6 x 10~^, we find reasonable agreement 
between the theoretical solution for £ and the experimen- 
tal data in Fig. [T5] 

We do see some disagreements between the experi- 
mental results in Fig. [151 especially at the 0.5 mM glu- 
cose concentration. In general, we expect to underesti- 
mate the penetration depth due to our approximation of 
colonies as hemispheres. The deviations may also indi- 
cate the presence of biological mechanisms that allow the 
cell clusters to adapt to particular nutrient environments. 
For example, yeast cells can modify their nutrient uptake 
parameter v = na by expressing different nutrient trans- 
porters on their surface [l|. Our theory assumes that 
each cell in a colony has the same absorption strength 
v, an approximation suitable provided all cells employ 
the same nutrient transporter system, and the ambient 
nutrient concentration is small enough so that we are 
in a linear regime of nutrient absorption in a Michaelis- 
Menten-type curve (see Fig. [5]). So, deviations from our 
theory may indicate deviations from the linear regime 
or the presence of regulatory mechanisms in cells that 
modify nutrient uptake kinetics during colony growth. 



CONCLUSION 
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Our theory also predicts that nutrient shielding is more 
sensitive to the volume fraction when each cell is a 
good absorber (i.e., for i/ > 0.1). As (p increases, cor- 
relations between the cells become more important as 
nutrient collisions with multiple cells create a stronger 
shielding effect. Thus, a colony of good absorbers, such 
as gram-negative bacteria, should be able to change its 
nutrient shielding properties by tuning the separation be- 
tween each cell. This hypothesis could be tested by plac- 
ing bacterial cells in a gelatinous matrix at various cell 
densities and observing their growth. As discussed in 
Sec. IIIDl these artificial colonies have many industrial 
and environmental applications. 

It would be interesting to extend our theory to include 
nonspherical cell shapes and cell radius polydispersivity. 
Polydispersivity should increase the screening length, but 
a detailed understanding of its effect on cell correlations 
is lacking. One could also account for spatial variability 
in (j). Cell clusters should "thin out" near the cluster 
surface, where the cells have not had time to grow into a 
densely packed structure. If this density variation occurs 
on scales large compared to the cell radii, it should be 
sufficient to replace </> by a spatially dependent 0(r) in 
our effective medium calculations. 

To extend the experimental results and test the effec- 
tive medium theory more precisely, one could control for 
the variability of the v parameter in the colony. One pos- 
sibility is to use engineered yeast strains with fixed nutri- 
ent uptake kinetics, such as the mutants constructed by 
Rcifcnberger ct al. [2[. These cells express a single type 
of glucose transporter in a medium with a low glucose 
concentration. In addition, it would be interesting to 
more precisely characterize the local glucose concentra- 
tion in the colony by using either a more direct reporter 
or more precisely characterising the relationship between 
glucose level and ribosomal protein expression. 



We have developed and tested an effective medium the- 
ory of nutrient transport in clusters of cells. Simulations 
and experiments support the theoretical results. The key 
parameters in our model are the single cell nutrient ab- 
sorption parameter v and the packing fraction </> of cells 
in the cluster. 

The wide variation in v for different cell types such as 
yeast (z/ ~ 10~^), gram-positive {f ~ 10~'^), and gram- 
negative {v ^ 10~^) bacteria has important implications 
for nutrient absorption. Effective medium theory pre- 
dicts very different nutrient shielding properties of cell 
clusters as i/ varies (see Figs. [TOl and [T2]) . Thus, the frac- 
tion of actively growing cells at a cluster surface will vary 
significantly with cell types. For instance, we predict that 
the screening length £ for an S. cerevisiae colony is about 
10 times longer than the length in a gram negative E. coli 
bacteria colony and 3 times longer than in a gram posi- 
tive L. mesenteroides bacteria colony, measured in units 
of the respective cell diameters. 
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Appendix A: Detailed Effective Medium Calculation 

Our effective medium calculation of the response func- 
tion E will closely follow Cukier and Freed's analysis of 
the perfectly absorbing cell case [42] . We generalize their 
work by allowing for partially absorbing cells via the ra- 
diation boundary condition (Eq. I13p . Cukier also consid- 
ered this case [43| , but only to first order in scattering in 
the effective medium. We extend his argument to include 
the pair distribution function of the cells, which leads to 
an important correction term. 

We begin by solving for the Lagrange multiplier func- 
tions crQ(ri) in Eq. [T2]in terms of an "effective medium 
propagator" Gw{j^,J^')- Upon assuming the effective nu- 
trient transport properties are homogeneous over the 
cell cluster, the propagator is translationally invariant 
(Gvy(r, r') = Giv(r — r')) and, in Fourier space. 



Gwici) = [Doq^ + W{q) 



(Al) 



where W is some suitable approximation to the self- 
energy E that tells us how nutrients diffuse through a 
medium consisting of the nutrient sinks and their inter- 
stitial space.. To simplify the analysis, we introduce an 
operator notation for all of our convolution integrals. For 
example, [W'^jr = J lF(r- r')^(r') dr'. Subtracting the 
convolution W^^ from both sides of the equation for ^(r) 
(Eq. [HI) yields 



*(r) 



dr'Gw{r-r') 



N 



-G 



w 






(A2) 



where CTi(r) = /dfJj (5(r — Rj)cri(ri), and Gh' [■ • ■] = 
/ Gw{y — r')[. . .]r' dr' is another example of the convo- 
lution described above. 

To solve for Ui via Eq. IA21 denote the inverse of the 
Gw operator on the i*^ ceh surface by K~^{ri, r^), where 
Ti and y[ are two vectors from the cell origin to the cell 
surface (see Fig.UJ^c)). Then, since Ea. lA2l is an equation 
for ^(r) for all points r in the cluster, ^(r) in Eq. IA2I 
is evaluated at an arbitrary point r = R^ on the i-th 
cell surface. In the perfectly absorbing case Cukier and 
Freed [42| considered, the boundary condition ^(R;) = 
is then used to solve Eq. IA2I for ai (r^ ) . The analogous 
condition for partially absorbing cells is Eq. [TS] Using 
the partially absorbing boundary condition yields 



<J,{r,) = -Jdn',g,{r,,r',) 



(i-q:)g 



w 



W* 




=R' 

(A3) 



where gi = [l — K- ^QiGw] Ki ^ is an operator de- 
fined on the surface of the i*^ sphere and does not de- 
pend explicitly on the sphere center R?. Q[ is the gra- 
dient operator with respect to the r' coordinate. Eq. IA3I 
corresponds to Eq. 2.9 in Cukier's analysis [i^. Note 
that in the operator notation, gi multiplying a func- 
tion implies a convolution over the i"^ sphere surface: 

g,[...] = /df}^g,(r„rO[...]r=R;. 

A useful "scattering" operator used in the following is 

r,(r, r') ^ I <5(r - R,)5.(r„ r',)(l - Q'J 

x5{v' -'R!i)dn,dn[. (A4) 

The operator Ti describes the scattering of a nutrient off 
of the surface of the i"^ cell. We now substitute Eq. IA3I 
into the di term in Eq. IA2I and solve Eq. IA2I for 'J by 
iteration: 



^ = Gw 



N 






w-^ 



i=l 



1 — Gy/T + y GwTiGy/Tj — . , 

X Gw{W^ -s) 

^ - [1 + GwiT -W)Y^ Gws, (A5) 



where T = J^i^^i- ^^ Eq. IA5[ it is important to in- 
corporate an exclusion condition for consecutive sums in 
the expansion of [1 + Gw(T - W)]-\ e.g. TGwT ^ 
Y.i j^i^i'^wTj (see ^^ for more details). Note also 
that Eq. [16] (main text) implies that the ensemble aver- 
aged field t/j = (^) in Fourier space satisfies 



- Doq^^ - Wt/J = Y.^^ - Wtp + s 
■ip^-Gw{Y-W)ip-Gws 

ip = -[l + Gw{Y- W)r^ Gws. 



(A6) 



The field -0 in Eq. IA6I must be equal to the ensemble 
average of Eq. IA5I We combine the two expressions and 
set W^ = E to find a self-consistent equation for E. After 
some algebraic manipulations (see Cukier's analysis [43| 
for details), the equation can be expanded in a series in 
T. The expansion up to second order is 

^^{T}- E (^'GsT,^,) + (T) Gs (T) . (A7) 

It makes sense to expand in T because this operator de- 
scribes a single interaction of the nutrient with any of 
the cells. The higher order terms in the expansion de- 
scribe multiple scattering events, which we expect to be 
less probable. We follow Cukier's analysis [43 in keeping 
only the first term on the right hand side of Eq. IA7I 

The various operators in Eg. I ATI arc computed by mov- 
ing to Fourier space and exploiting expansions in spher- 
ical harmonics. For example, Gw{^i,^i) = Gwi^i — 
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rj) (given by Eq. lAip is expanded in Fourier modes 
giq(ri-r.)^ which Ere then rewritten in terms of spher- 
ical harmonics via the spherical wave expansion of the 
plane wave [S^- We find 

Gw{^^,^';) =Y^iiY,mmYl„^{^,), (A8) 



where Z = 0, 1, . 



m = —£, —i+l,...,£, and 

2q^dq iiiqaf 
IT Doq^ + Wiq) ■ 



(A9) 



The ji functions are the spherical Bessel functions of the 
first kind. Cukier [43| describes this procedure in more 
detail. Another important result (Cukier's Eq. A15) is 

i,m 

(AlO) 
where 



. ^ 2a f^ q^f,{qa)Mqa) 
"-' - ^v Jo '^ Doq^ + Wiq) 



(All) 



Now that we have an expression for g,;, it is possible to 
compute [T). The bracket ensemble average will require 
an averaging over all cell positions, so let us assume that a 
single cell center is distributed uniformly over the cluster 
volume V, so that P(R°) = V^^ for alH = 1, . . . ,iV. 
After averaging, the operator (T) acts on an arbitrary 
function /(q) in Fourier space as follows: 

(T) /(q) = c f e-'^-^o^oCro, r^)(l - ^^-'r', ■ q') 



X e 



'q'-ropiCq'-qj-r 



/(qOdf^od^ 



, dq' 



(27r)3 



dr 



= 47rc^K£(2£+l)jf(ag) 

e 

X [jeim) - J^^^aqfiiai)] /(q) 

^r(q,q)/(q). 



(A12) 



where c = N/V is the concentration of cells and 
Kg = [1 - 77^0]"W7^ In the last line of Eq. \M% 
we have introduced T(q, q'), a convenient scattering op- 
erator (averaged over the single cell center distribution 
P(R^) = V^^). In Fourier space, 

7"(q, q') EE 47r^ K£ [Xag') - v^^ aci j[{a(i)\ 
f. 

xji{aq){2e + l)Pi{cose,), (A13) 

where 9q is the angle between q and q' and Pa are the Leg- 
endre polynomials. Only the diagonal elements T(q, q) 
contribute to <^T). The off-diagonal terms will contribute 
to the second order correction term ^^ -,^ (TiG^Tj^i) in 
Eq.E3 



The second order correction term requires the use of 
the pair distribution function P(R°,Rp for the centers 
of pairs of cells i and j. We use the approximation Ea.[3T] 
for P(R0,R0) and find 
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X /r(q,q')GE(q')^(q'-q)r(q',q) '^'^' 



(27r)3' 
(A14) 



where ^(q) is the Fourier transform of 6'(|r| — 2a): 

^ 32^a3M2a^^ 

Zaq 



(A15) 



We assume N is large enough so that N{N — 1) w A^^, 
and then substitute Eq. IA14I into Eq. IA7I to find that to 
second order in scattering in the effective medium 

S(q)-cr(q,q) 

4a-'^c2 f . ,^.7i(2a|q'-q|) , , 

+ ^^ y nq, q )G^(q ) ^^^^^:^r(q , q) dq . 

(A16) 



In principle, Eq. IA16I could be solved for S(q) numer- 
ically using an iterative procedure. However, the sec- 
ond order terms involve an integral of G^{q) = {Doq^ + 
T,{q))^^ multiplied by various Bessel functions. Due to 
the oscillatory nature of the Bessel functions, these kinds 
of integrals have poor convergence properties using stan- 
dard numerical techniques and require special integration 
methods [57[. To avoid these cornplications, we make 
a "hydrodynamic" approximation [32, 43, 58] in which 
I](q) is approximated by the first two terms of its Taylor 
expansion: S(q) ^ k + SDq^. 

The constants 7^ and Q have closed forms in the hy- 
drodynamic approximation and involve integrations over 
spherical Bessel functions (tabulated in [53 )• We first 
note that 



7f 



Dna 



kixf 



■ X dx ~ 



aii{a)ki{a) 



(A17) 



where a = a/^ is the ratio of the cell radius to the corre- 
lation length ^ = y^D/k and ii{x), ki{x) are the modified 
spherical Bessel functions of the first and second kinds, 
respectively. Next, we have 



TTvDa 



■M^±M^x'dx 



a'^i£{a)k'g{a) 



uDa 



(A18) 



where the x^ means wc evaluate the derivative of the 
Bessel function at x-l-e and let e ^- 0. Using Ea. lA17l and 
Ea. lA18l to evaluate T(q, q') in Eq. lA13l and substituting 
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T(q, q') into Eq. IA16I (evaluated at q = 0) yields 



cr(0, 0) 



8aV 
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ATTDaai'{l + cotha) 



1 2 
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jo{x)ji{2x) [jo{x) - V ^a;Jo(x)] 



sdx, 
(A19) 



where the integral over the product of three spherical 
Bessel functions is performed by writing the Bessel func- 
tions in terms of elementary functions such as powers and 
exponentials and then using contour integration tech- 
niques. The final result reduces to the implicit equation 
for ^, given in the main text as Eq. 1321 

The diffusion coefficient D can be calculated by look- 
ing at the second order term in the Taylor expansion 
of Eq. IA16I around q = 0. Note that for the purposes 
of determining the effects on nutrient shielding in a cell 
cluster, the relevant quantity is the screening length ^. 
A lengthy calculation of the change in the diffusion coef- 
ficient bD = D — Dq yields 
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where a = a/£, is given by the self-consistent solution 



to Eq. [32] and [ke{a)Y 
£=0,1. 
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